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Inelastic collisions occur in Bose-Einstein condensates, in some cases, producing particle loss in the 
system. Nevertheless, these processes have not been studied in the case when particles do not escape 
the trap. We show that such inelastic processes are relevant in quantum properties of the system such 
as the evolution of the relative population, the self trapping effect and the probability distribution 
of particles. Moreover, including inelastic terms in the model of the two-mode condensate allows 
for an exact analytical solution. Using this solution, we show that collisions favor the generation of 
entanglement between the modes of the condensate as long as the collision rate does not exceed the 
natural frequency of the system. 



I. INTRODUCTION 

Bose-Einstein condensates are quantum systems of 
great interest since they consist of up to 10 10 particles 
which exhibit collective quantum behavior. Recently, 
multi-component (or multi-mode) Bose-Einstein conden- 
sates have been extensively studied both theoretically 
QL S Hi and experimentally 0, [ESQ- The components 
of the condensate are either several spatially separated 
condensates Q or several internal degrees of freedom of 
the condensed particles 0). Unfortunately, the lack of 
analytical solutions for these systems limit our under- 
standing of them. Until now, multi-component conden- 
sates have only been studied through numerical solutions 
which are always limited by the growing degrees of free- 
dom of the system. 

The simplest multi-component condensate is known 
as the two-mode Bose-Einstein condensate. The sys- 
tem consists of condensed particles in two hyperfine lev- 
els [l|. Transitions between the levels are induced via 
a Josephson-type interaction produced by a laser. Two 
body collisions are, together with the laser interaction, 
the most relevant physical process in the system. An- 
other incarnation of the two-mode condensate consists 
of two spatially separated condensates coupled through 
a tunneling barrier corresponding to the Josepshon-type 
interaction (j|. The canonical model used to describe 
these systems considers the Josephson-type interaction 
and elastic collisions between same and different particle 
type [T[ . The ground state of the model can be found in 
the first or second quantization only by numerical means 

In a recent publication, we introduced a model for the 
two-mode Bose-Einstein condensate which includes in- 
elastic collisions 9]. In our model the inelastic collisions 
are induced by the Josephson type-interaction. Including 
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these terms allows for an exact analytical solution of the 
model. Several types of inelastic collisions are known to 
be relevant in real physical situations [T(J El, [HJ ■ The 
inelastic process better known in experiments are those 
which result in particle loss such as background collisions, 
three-body recombination, spin-exchange and dipolar re- 
laxation [TT|, [l3| . The inelastic processes that we consider 
here have not been properly characterized by experiment 
yet. The main reason for this is that they do not give 
rise to particle loss. Nevertheless, there is good indica- 
tion that this process are present in real situations. In 
this paper, we present a number of theoretical predic- 
tions on the evolution of the relative population of the 
condensate, the self trapping effect and on the probabil- 
ity distribution of particles in the condensates to be cor- 
roborated by experiment. We compare our results with 
the predictions of the canonical two mode Bose-Einstein 
condensate. Experiments could be performed to compare 
the predictions of the models. The great advantage of our 
model is that it has exact analytical solution which allows 
for a complete understanding of the system. Remarkably, 
we find that the ground state of the system is the coher- 
ent state which was found to provide a good description 
of the two-mode Bose-Einstein condensate [IJ, [l5| . We 
take advantage of the simplicity of the model to show un- 
der what conditions maximally entangled states between 
the modes can be created and the role of collisions in the 
generation of entanglement. 



II. PARTICLE LOSS DUE TO INELASTIC 
COLLISIONS IN BOSE-EINSTEIN 
CONDENSATES 

Inelastic collision in Bose-Eistein condensates occur 
with high frequency and have been studied extensively 
[ToL [ill O, [TH . These processes are well known since they 
commonly result in particle loss. For this reason there 
has been considerable efforts to suppress them in exper- 
imental situations (Til Il3| . The main inelastic processes 
known are background collisions, three-body recombina- 
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tion, spin exchange and dipole relaxation. Atoms in a 
condensate are commonly knocked out of the trap by 
collisions with untrapped, room-temperature molecules 
in the background chamber. In the three body recombi- 
nation process three particles collide and form a molecule 
which is no longer trapped by the potential. Particles are 
also lost when, during a collision, the atoms make transi- 
tions to different energy levels which are not trapped by 
the potential. Spin-exchange and dipolar relaxation are 
important mechanism of particle loss which are closely 
related to the process considered in our model. After a 
collision either or both atoms exit in a different spin state. 
If the spin re-orientation energy is larger than the depth 
of the trapping potential this mechanism might result in 
the loss of the particles from the trap [TB] ■ In dipolar re- 
laxation one or two particles in the higher hyperfine level 
decay to the lower level during the collision due to their 
dipole-dipole interaction. The excess of energy is then 
transferred to an excess of momentum which produces 
the particle to escape the trap. In our model the inelas- 
tic process occurs but the excess of energy is not enough 
to help the particle escape from the trap. Theoretically, 
inelastic processes are included in the studies of the con- 
densates by considering classical rate equations for the 
particle loss . But in the literature there is a lack of 
experimental and theoretical analysis on inelastic process 
which do not give rise to particle loss. Our model exem- 
plifies the importance that these processes might have in 
the behavior of the system. 



III. TWO-MODE BOSE-EINSTIEN 
CONDENSATE WITH INELASTIC COLLISIONS 

Our model of the two mode Bose-Einstein condensate 
consists of the following Hamiltonian, 

H 2 = A +5w((Ja-tfb) (1) 

+ \{e^a^b + e- ilt >atf)+U^tfab 

+ A(e 2i Va^fe + /i.e.) 

+ Li{{a)a)ab-tfa)bb)e i ' t ' + h.c.). 

The modes a\ a and b\ b with frequency difference 5ui, 
correspond to either atoms with two different hyperfine 
levels or alternatively, two spatially separated con- 
densates [H . The Josephson-type interactions is induced 
by applying a laser [5| or a magnetic field gradient 
In our Hamiltonian the Josephson-type term, in which 
one particle is annihilated in one mode and created in 
the other, has coupling constant A and phase (p. The 
terms with four bosonic operators describe two-particle 
clastic and inelastic collisions. The elastic collisions have 
interaction strength U. The inelastic collisions have in- 
teraction strength fi when two particles in the same mode 
collide and one of them is transformed into the other and 
interaction strength A when the collision transforms two 



particles in one mode into the other. Ignoring the in- 
elastic terms in Eq.(TT]), we find that our model coincides 
with the canonical Josephson Hamiltonian [l[ when the 
rate of collisions of same particle type is equal. The as- 
sumption of equal collision rates for the same particle 
type is also made in [l[ in order to find approximate and 
numerical solutions. Our model has the same number of 
free parameters as the canonical two-mode Hamiltonian. 
The only difference is that the Hamiltonian in Eq. (JTJ) in- 
cludes inelastic collisions, which are usually present in 
realBECs 

Including the correct rate of inelastic collisions in the 
Hamiltonian allows for an analytical solution. The solu- 
tion is simply found by considering the following Hamil- 
tonian Hq — Ai((ra — b^b) + A 2 (a)a — b^b) 2 with eigen- 
states \N, m) where N is the eigenvalue of the total num- 
ber operator N = n a + rib = (V a + wb and m the eigen- 
value of the relative population rh — a*a — b'b. Since the 
number of particles in the system N is constant, m is 
restricted to values m = — N, ...,N. Ho describes a two- 
mode condensate with no Josephson-type interaction, i.c, 
no laser interaction or tunneling barrier. The parameter 
A 2 corresponds the two-body elastic scattering probabil- 
ity and A\ to the energy difference between the modes. 
It is easy to verify that if the parameters of the model 
satisfy the following conditions 



A = 


A 2 (cos 2 9N 2 +sin 2 ON)/ A, 


Slo = 


{A x cos6>)/2, 


A = 


(Ax sin0)/2 


U = 


A 2 {\ -3 cos 2 0)/4, 


n = 


(A 2 cos6»sin6>)/2, 


A = 


(A 2 sin 2 6) /A. 



then the Hamiltonian can be wr itten as H 2 = UH W 
where where U is the two-mode displacement operator 
U = e^ b+ t afct with displacement parameter £ = 9e 1 ^. 
The solution of the Hamiltonian (JTJ is simply W\N, m) 
with energy £ m = Axm 4- A 2 m 2 . 

Note that by fixing the free parameters 5u>, A and U 
the inelastic collision constants, A and fi, are determined. 
This is because in our model, the inelastic collisions are 
produced by the effect of the Josephson-type interaction 
on colliding particles. This physical relationship is math- 
ematically expressed by the relationship between the co- 
efficients. The relationship between the parameters im- 
ply that, for fixed 9, the scattering probabilities for in- 
elastic collisions are higher if the scattering for elastic 
collisions are more probable. If sin(0) = 0, the Joseph- 
son coupling is zero and there is no inelastic processes 
occurring in the system. There are only elastic collisions 
between particles of the same type. We also see that 
for a resonant laser or a symmetric well 5lo = 0, which 
corresponds to cos 9 = 0, there are elastic and inelas- 
tic collisions but the inelastic interaction always involves 
the exchange of two particles simultaneously. We would 
like to emphasize that if for a specific physical system, 
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the rates of inelastic to elastic collisions does not hold or 
cannot be arranged by external manipulation, an analyt- 
ical solution cannot be found using our method for such 
a condensate. Fortunately, in the laboratory the rate of 
elastic and inelastic collisions can be manipulated, for 
example, by applying a magnetic potential [171 ] making 
possible to meet the experimental values for the produc- 
tion ratios of those terms. The advantage of having a 
region in the parameter space of the system where we 
can find an exact solution is that it is then possible to 
find solutions for other parameters using perturbation 
theory. 

Due to the simplicity of our solution the ground state 
U'\N, mo) of H 2 is trivially found by minimizing the 
energy — Aim + A 2 m 2 with respect to m. For 
A 2 > 0, rrio is the nearest integer to —Ai/(2A 2 ) or 
to = -AiiV/|Ai| when | -A 1 /(2A 2 )\ > N. For A 2 < 0, 
the minimum corresponds to too = N if A% < or 
to = — N otherwise. The solution for A 2 < and Ai > 
is the spin coherent state which has been argued to de- 
scribe the Bose-Einstein condensate well enough [l4|, EH • 

IV. INELASTIC COLLISIONS IN A DOUBLE 
WELL 

In the physical realization of the model, where the two 
modes a, and b, b^ of the condensate correspond to par- 
ticles in two spatially separated condensates, A and B, 
the Josephson-type interaction corresponds to the tun- 
neling of particles through a barrier with probability pro- 
portional to A. 

The canonical model, as well as our model, considers 
elastic collisions between the particles in each well given 
by the terms proportional to a^aa^a and Two 
particles in one well are annihilated and two particles in 
the same well are created. Note that in our model these 
terms are part of the constant term N which is included 
in the parameter Aq. Both models also consider colli- 
sions between particles in different wells given by a^ab^b. 
One particle in well A and one particle in well B are an- 
nihilated and one particle in well A and one particle in 
well B are created. These collisions occur in the region 
where the spatial wave functions of each well overlap. In 
all these elastic interactions the particle number in each 
well is conserved. The new feature in our model is that 
we consider collisions occurring in the overlapping region 
in which particles in a well end up in the other well after 
the collision. We consider terms in which two particles 
in one well collide (close to the tunneling barrier) and 
one or two of them are simultaneously tunneled to the 
second well. We also include situations in which parti- 
cles from different wells collide in the overlapping region 
and both of them end up in the same well after the colli- 
sion. The probability of such inelastic events depends on 
the overlapping region between the condensates, which 



is a function of the tunneling probability and the trap- 
ping potential. The probability of inelastic events is also 
a function of the probability of elastic events which de- 
pend on the scattering lengths of particles in the system. 
It seams logical and natural to consider such inelastic 
processes. Note that there are other theoretical motiva- 
tions to support the inclusion of such terms 18]. The 
specific relationship between the parameters for which 
the model has an exact analytical solution corresponds 
to a special type of potential. 



V. INELASTIC COLLISIONS BETWEEN 
PARTICLES IN THE PRESENCE OF A LASER 
FIELD 

In the physical realization of the model where the two 
modes correspond to two hyperfine levels, the Josephson- 
type interaction corresponds to the interaction of a laser 
field which induces transitions between the atoms hy- 
perfine levels. In these picture, same and different spice 
particle can collide. The inelastic collisions correspond to 
the collision of particles in the presence of a laser field. 
Energy is absorbed and liberated into the field during col- 
lision in which the particles change their hyperfine levels. 
Although a detailed microscopic calculation of this pro- 
cess has not been yet calculated, it has been observed in 
experiments that the rate of inelastic collisions between 
atoms is increased when there is an interaction with a 
laser field [\ '1 . We provide in this paper several theoret- 
ical predictions to be experimentally tested of the effects 
of inelastic collusion which do not give rise to particle 
loss in the condensate. 



VI. EFFECTS OF INELASTIC COLLISIONS IN 
THE SYSTEM 

In this section, we analyze the effects of the inelastic 
collision terms added to the Hamiltonian in the behavior 
of the two-mode Bose-Einstein condensate. We study the 
particle distribution of the ground state and the evolu- 
tion of relative population in each mode. We will show 
that the presence of inelastic collisions have important 
consequences in the self trapping effect. 



A. Particle distribution of the ground state 

Let us first consider the analytical solution to our 
model. The particle distribution P = \(N, m^o)! 2 of 
the ground state | V'o) = U'\N,mo) is equal to the mod- 
ulus square of the Clebsch-Gordon coefficients (Wigner) 
P=\d%, mo \ 2 where, 
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i N 



y/(N + m )\(N - m )\(N + m)l(N - to)! 
(N + m - k)lkl(N — k — m)\(k - m + m)\ 



(cos(0/2)) 



2N-2k+m -m 



(sin(0/2)) 



(3) 



The sum must be done over k whenever none of the argu- 
ments of factorials in the denominator are negative. Dif- 
ferent ground states parameterized by Too are obtained 
by changing the rate Ai/A 2 since to is the closest integer 
to -A l /2A 2 (for | - A 1 /2A 2 \ < N otherwise, m Q = N). 
When toq — N, the particle distribution is 



2NI 



(N + m)\(N - to) 



cos(6y2) JV+m sin((9/2) 



N—m 



which corresponds to a distribution with one maximum. 
We plot a n example for N = 1000 particles and = 1 in 
Fig. 1(a) As we change the rate Ai/2A 2 from 1000 to 
977, the distribution changes abruptly from a distribu- 
tion with a single maximum to a distribution with sev- 
eral peaks corresponding to multiple cat states as shown 
in Fig. [TJ These discontinuous changes occur when the 
system undergoes first order transitions corresponding to 
level crossings. The structure of the ground state changes 
abruptly and this is reflected in the probability distribu- 
tion. As too approaches zero, the number of maxima in 
the distributions grows. Controlling the inelastic param- 
eters allows us to prepare the condensate in such super- 
position states. 
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distribution with one maxima to a superposition distri- 
bution is no longer abrupt. We find the abrupt transi- 
tion only when the parameters change along the values 
corresponding to those for which the system has exact 
analytical solution. In the canonical model, which does 
not include inelastic collisions, the change from a sin- 
gle peaked distribution to a double superposition is also 
smooth as shown in [l| . The main difference observed be- 
tween the canonical model and our model is that in the 
canonical model only single and double superposition dis- 
tributions are observed. In our model, it is possible to 
obtain the superposition of N distributions. Therefore, 
the presence of more than two components in the prob- 
ability distribution is a clear footprint of the presence of 
inelastic collisions in the system. 
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FIG. 2: Ground state relative population distribution for 1000 
atoms. 5uj = 109 J, A = 487J and U = 0.214027J. 



FIG. 1: Ground state relative population distribution for 1000 
atoms and 9 = 1. Different mo correspond to a different rate 
between the parameters Ai and A2. Quantum superposition 
appears when mo < 1000. 



B. Inelastic collisions in the evolution of the 
relative population and the self trapping effect 



Now consider the Hamiltonian ([I]) with general pa- 
rameters. In this general case it is not possible to find 
an analytical solution. Nevertheless we find, by numeri- 
cal means, the probability distribution of particles in the 
ground state. In Fig. [2] we show the probability distribu- 
tion for N = 1000. We leave 5u> = 109J, A = 487J and 
U = 0.214027J fixed and vary \i and A as shown in the 
figure. By varying these parameters we show that the 
transition from different probability distributions of the 
ground state can be done smoothly. The change from a 



We now calculate the evolution of the average relative 
population (a^a — b^b), for the initial condition \ip(t = 0)). 
In the case where the system has an analytical solution, 
the evolution of the relative population for a given initial 
state \tp(t = 0)) = T,m=-N C mU^\N,m), with coeffi- 
cients C m , is given by 

JV 

(a f a - tfb) = cos 0^2 to \C m | 2 (5) 

— JV 
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N 

- sin6 C m C m -x(N(N +l)-m(m-l)) 1 / 2 L m 

-N+l 

L m = cos(0 + (-E m _i - £ m ) t) (6) 

In the general case, the evolution of the relative pop- 
ulation is calculated numerically. In Fig. ([3]) we plot 
the evolution of the relative population for Suj = 109 J, 
A = 487 J and U = 0.2 14027 J and vary fj, and A as 
shown in the figure. Fig 3 (a) correspond to the evolu- 



tion of the canonical model (without inelastic collisions). 
In Figs.( 3(b)]|3(c) I we plot the evolution of the general 
model which includes inelastic terms. We considered sev- 
eral cases with different amount of inelastic collisions. 
Until here all the results where numerical. Finally, the 
exact analytical solution given by Eq. (J5j) is plotted in 
Fig |3(d)| In the canonical model the collapse and re- 
vivals are not uniform. Each revival becomes wider and 
at some point one revival gets mixed up with the next 
one. We can see that adding inelastic collisions changes 
drastically the behavior of the system. As the rate of 
elastic and inelastic collision approaches the values for 
which we have analytical solutions, we observe that the 
revivals start becoming more uniform. We can go from 
the solution of the canonical model to the analytical so- 
lution of our model continuously by varying the inelastic 
collision parameters. Finally, for the analytical solution 
of the model, obtained when the rate of elastic to inelas- 
tic collisions is equal to 1/2, the evolution of the rela- 
tive population is periodic and all revivals have the same 
width. One can see in Fig |3(d)| that the structure of the 
revivals repeats after two collapses. 

We take advantage of having an analytical solution to 
calculate (using Eq. ([5])) the time at which the oscilla- 
tions of the relative population collapse and the period 
at which the collapse and revivals occur. In order to cal- 
culate the time of collapse we find the condition when 
the different terms in Eq. ([5]) cancel each other. This 
occurs when the phase difference between the cosines in 
L m is (2n + l)n, where n is an integer. Therefore, we 
obtain the condition (-E m _i — E m )t — (E m — E m+ i)t = 
7r(2n -(- 1) which implies that the time of collapse is 
t r = (2n + l)ir / (2A2) . We will assume that t = cor- 
responds to exactly half of the revival (as in Fig. [3]). 
Therefore, t r (n = 0) corresponds to half the distance be- 
tween the revivals. 

To calculate the period of the evolution of the relative 
population we suppose that the initial time is to = and 
consider a later time t\ > to &t which the average relative 
population is the same {a) a — b%){to) = (a^a — a^a)(ti). 
If for every integer m exists an integer n m where the 
following condition is satisfied, 



(E m -i - E m )t x = 2-K<n r , 



(7) 



from Eq. ([5]) we can see that the evolution of all the terms 
L m is identical from t\ as it is from to- If we consider 
t\ to be the smallest number which satisfies Eq. ([7]) the 
function (a^a — b^b) is periodic with period t\. We can 
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FIG. 3: Evolution of the relative population (m),5uj = 109 J, 
A = 487J and U — 0.214027J and vary fj, and A are shown in 
the figure. 



write Eq. ([7]) as, 

Mi 



A 2 (2m-l))ti 



27T7V, 



(8) 



If A1/A2 is not a rational number Eq. (jHJ) is never sat- 
isfied and the function is not periodic. We will have 
collapse and revivals, but the particular form of every 
revival would be different. When A1/A2 is a ratio- 
nal number, the time t\ corresponds to the period, and 
pr = t\/{2t r {n = 0)) is the period of the relative pop- 
ulation in units of number of revivals. Remember that 
t r (n = 0) corresponds to half the distances between re- 
vivals. In other words, p r — k means that after k revivals 
the particular structure of the first revival repeats itself. 
Considering A\ — cA^ with c = p/q and p,q integers, 
one can se that if —q+p is even then pr = q, and if 
—q + p is odd then pr — 2q. An example is shown in 
FiglU where the evolution of the relative population is 
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plotted for the following values 9 = 1.35, A2 — 1 and 
Al = 49,50, 101/3,59/2. The values are chosen so that 
Pr = 1,2,3,4. 

In the canonical model (jj, = 0, A = 0), there are no 
Rabi oscillations if W > A. The condensate remains 
trapped in one mode if the scattering probability is much 
larger than the Josephson-type interaction. This corre- 
sponds to the well-known self trapping effect [lj|. Now 
consider our model, when the parameters satisfy condi- 
tions ((2]) and the evolution of the relative population is 
given by Eq. ([5]). If A 2 and Ai are fixed, the ratio U/\ is 
maximum for 8 = 0. Under these circumstances, it can be 
trivially seen from Eq. ([5]) that there are no oscillations 
and thus, we observe self trapping. This case is trivial 
since there is no Josephson-type interaction. However, if 
we consider 9 fixed, we find that the inelastic collisions 
can prevent self trapping from occurring even if Li 3> A. 
For instance, consider 9 = tt/2, there is a collapse at 
time t r (defined above). In this case, although A 2 3> A\ 
(implying that U ^> A) there is never self trapping. 



VII. ENTANGLEMENT 

Entanglement is a property of multipartite quantum 
states that arises from the tensor product structure of the 
Hilbert space and the superposition principle. It plays an 
important role in the understanding of many body quan- 
tum systems. Moreover, it is considered to be a resource 
for quantum information tasks [201 ] such as teleportation 
and has applications in quantum control and quantum 
simulations [2l|. It is therefore of interest to find the 
necessary conditions to prepare states with high degrees 
entanglement. 

Entanglement between the modes of a two component 
Bose-Einstein condensate have been studied in the canon- 
ical model [TBI, [22j • In this section, we calculate entan- 
glement between the modes using the analytical solution 
to our model and show the conditions to generate max- 
imal entanglement in the system. We analyze in detail 
the role of collisions in the generation of entanglement. 

Quantifying entanglement for pure bipartite systems 
is simple. Entanglement between the systems is given 
by the von-Neumann entropy S(p a ) = — tr(p a log 2 (p a )), 
where p a is the reduced density matrix. The entangle- 
ment between the modes for the ground state of our 
system W\N,mo) is given simply by S(9, N,m ) = 
-EMd% ,„M\ 2 log 2 \d? n(um (9)\ 2 where d^ m (9) are 
the Clebsch- Gordon coefficients. The entanglement is a 
function of the angle 9, the number of particles N and 
mo. The ground state is determined by the parameter 
mo which depends on the rate Ai/2A 2 . This parameter 
regulates the size of the elastic collision probability in 
comparison to the difference of frequency of the modes 
in the absence of Josephson-type interaction. 

It is easy to see from the formula that the entangle- 
ment of the state mo = N is maximum when 9 = tt/2 
independently of the number of particles in the conden- 



sate. This means that when the collisions rate is much 
smaller than the frequency difference between the modes 
in the absence of Josephson-type interaction A 2 -C A\, 
the entanglement generated after turning on the interac- 
tion, will be maximum if the double well is symmetric 
with a high tunneling probability. In a two species Bose- 
Einstein condensate, maximum degrees of entanglement 
will be generated, in this case, by turning on an intense 
resonant laser. Note that when mo = -/V the probability 
distribution has a single component. 

We plot in Fig.fJS^,) the entanglement between the 
modes for different ground states mo as a function of 
9 for N = 1000 particles. We find that the entanglement 
is generally higher when m is small. This corresponds to 
a regime where the probability distribution of particles 
has many maximum peaks. This occurs when the elastic 
collision rate becomes comparable to the frequency dif- 
ference between the modes in the absence of Josephson- 
type interaction, A\ ~ A 2 . But note that by turning on 
the Josephson interaction with 9 close to tt/2, there is 
a local minimum at mo = 0. In Fig. (JSJa) we plot en- 
tanglement at 9 — tt/2 and observe this local minimum. 
This means that for a large Josephson-type interaction 
strength, there is a local minimum when A 2 > A±. 

We have learned something interesting about the role 
of two-body collisions in the generation of entanglement 
in a Bose-Einstein condensate. Higher degrees of entan- 
glement are generated by switching the Josephson-type 
interaction when the rate of collisions is comparable to 
the frequency difference between the modes (in the ab- 
sence of the Josepson-type interaction) A\ « A 2 . How- 
ever, if A 2 > A±, lower degrees of entanglement are 
generated. In the case of a double well condensate, if 
the collision probability in the absence of tunneling is 
comparable (but lower) to the energy difference between 
the wells, then, by letting the atoms tunnel lowering the 
barrier maximally, entanglement will be maximum if the 
wells are symmetric. However, if the collision probability 
exceeds the energy difference between the wells in the ab- 
sence of tunneling, maximum degrees of entanglement are 
achieved when lowering maximally the barrier if the wells 
are slightly asymmetric. This will favor the generation 
of entanglement. In the case of a Bose-Einstein conden- 
sate consisting of two internal degrees of freedom, if col- 
lisions occur with probability comparable (but lower) to 
the energy difference between the modes, entanglement 
will become maximum by coupling the internal degrees 
with a strong resonant laser. But if the collision proba- 
bility exceeds the energy difference between the internal 
levels (in the absence of the laser), it becomes more con- 
venient for the generation of entanglement, to turn on 
a slightly detuned laser coupling. The same conclusions 
can we drawn from the analysis of Fig. (|5^) where we 
plot the entanglement at mo = 0, i.e. A% > Ai. We ob- 
serve local minimum at 9 = tt/2. Therefore, we conclude 
that in this case, it is more convenient to have 9 sightly 
different to tt/2. 

In Fig. (T5Jd) we plot the entanglement as a function of 
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FIG. 4: The evolution of the relative population (m) is plotted for the following values theta = 1.35, A2 = 1 and Al 
49, 50, 101/3, 59/2. The values are chosen so p r = 1, 2, 3, 4 



the number of particles TV and different ground states mo 
for (9 = 7r/2. We can see that the entanglement is stronger 
as we have more particles in the condensate. We plotted 
the entanglement as a function of the number of particles 
for N as big as 10 5 and found the same behavior. 

Therefore, to create high degrees of entanglement be- 
tween the modes, it is convenient to have a condensate 
with large number of particles where the two-body elas- 
tic collision rate is comparable (but not larger) to the 



frequency difference between the modes in the absence 
of the Josephson-type interaction. One then turns on a 
strong near resonant laser (or lower maximally the tun- 
neling barrier in an almost symmetric well). 
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FIG. 5: Entanglement a) as a funtion of 6 and mo for N = 100 
and b) as a funtion of N and mo = for = tt/2. 
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FIG. 6: Entanglement for N=100 as a function of a) 9 for 
mo = and b) mo for 9 = tt/2. 

VIII. DISCUSSION AND OUTLOOK 

We have presented several predictions of a new model 
describing a two-mode Bose-Einstein condensate which 
includes spin-exchange inelastic collisions. These colli- 
sions have been observed in experiments when they lead 
to particle loss. Nevertheless, there are good reasons to 
believe that these processes do not always lead to particle 
loss. The particles escape the trap after an inelastic pro- 
cess only when there is an excess of energy larger than the 
trapping potential. Our work is the first to analyze the 



effects of inelastic collisions in the system when the parti- 
cles do not escape the trap. We show that such inelastic 
collisions produce important qualitative and quantitative 
effects. Moreover, our model has the great advantage of 
having exact analytical solutions for a set of parameters. 
Previously, the study of the two-mode Bose-Einstein con- 
densate was restricted to numerical treatment. 

Inelastic collisions which do not lead to particle loss 
have not been yet a topic of experimental research. The 
main intension of our work is to present a number of pre- 
dictions in the evolution of the relative population of the 
condensate and in the probability distribution of particles 
in the system which would allow to detect such inelastic 
process in the laboratory. The presence of these effects 
can be easily detected by comparing the results to the 
canonical two-mode hamiltonian which does not include 
inelastic collisions. We find that the presence of inelastic 
collisions produce multicomponent superposition states, 
rather than two component cat states as predicted by the 
canonical model. We show how inelastic processes have 
qualitative effects in evolution of the relative population 
and are capable of preventing the self trapping effect from 
occurring, even when the collision rate is larger than the 
Josephson-type coupling. Taking advantage of the ana- 
lytical solution of the system, we include in our work a 
study on mode entanglement and discuss under what con- 
ditions high degrees of entanglement can be generated. 
We find that higher degrees of entanglement can be gen- 
erated when the condensate consists of large number of 
particles and the two-body collision probability is com- 
parable (but not larger) to the energy difference between 
the modes in absence of Josephson-type interaction. In 
such condensates, maximally degrees of entanglement are 
generated by a strong resonant laser coupling (or small 
symmetric tunneling barrier in the case of a double well 
condensate). The analytical solution to our model helped 
us learn details about the role of two-body collisions in 
the generation of entanglement in the system. We found 
that collisions favor entanglement as long as the colli- 
sion probability does not exceed the frequency difference 
between the modes (in the absence of Josephson-type in- 
teraction). If such is the case, higher degrees of entan- 
glement can achieved by allowing maximal tunneling in 
a slightly asymmetric well (or turning on a strong but 
slightly detuned laser coupling). 

To complete our study of inelastic spin-exchange colli- 
sions in the two-mode Bose-Einstein condensate, we plan 
to introduce the loss of particles due to inelastic collisions 
as a decoherence process using a master equation. We are 
currently working in the generalization of our model to 
spin- 1 (or three well) condensates. We are also studding 
the effects of many-body collisions in the system. 

We thank C. G. Hernandez-Salinas for his help verify- 
ing some of our results. 
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